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Abstract: The simplest supersymmetric model that solves the fi problem and in which 
the GUT-scale parameters need not be finely tuned in order to predict the correct value of 
the Z boson mass at low scales is the Next-to-Minimal Supersymmetric Standard Model 
(NMSSM). However, in order that fine tuning be absent, the lightest CP-even Higgs boson 
h should have mass ~ 100 GeV and SM couplings to gauge bosons and fermions. The only 
way that this can be consistent with LEP limits is if h decays primarily via h — > aa — > 
t + t~t + t~ or 4j but not 46, where a is the lighter of the two pseudo-scalar Higgses that are 
present in the NMSSM. Interestingly, m a < is natural in the NMSSM with m a > 2m r 
somewhat preferred. Thus, h — > t + t~t + t~ becomes a key mode of interest. Meanwhile, 
all other Higgs bosons of the NMSSM are typically quite heavy. Detection of any of the 
NMSSM Higgs bosons at the LHC in this preferred scenario will be very challenging using 
conventional channels. In this paper, we demonstrate that the h — > aa — > t + t~t + t~ decay 
chain should be visible if the Higgs is produced in the process pp — > p + h + p with the 
final state protons being measured using suitably installed forward detectors. Moreover, 
we show that the mass of both the h and the a can be determined on an event-by-event 
basis. 
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1. Introduction 

The Next-to-Minimal Supersymmetric Standard Model (NMSSM) extends the MSSM by 
the introduction of a singlet superfield, S. To do so provides the possibility to solve [1] the 
fine tuning problem that arises in the case of the MSSM and it provides a natural solution 
to the ^-problem via the XSH U H^ superpotential term when the scalar component of 5 
acquires a vev ([J, e ff = A(S')). The Higgs sector of the NMSSM contains three CP-even 
and two CP-odd neutral Higgs bosons, and a charged Higgs boson. The lightest CP-even 
Higgs boson is denoted by h and the lightest CP-odd Higgs boson is denoted by a. 

In [2] and [3], a partial 'no-lose' theorem for NMSSM Higgs discovery was established. 
The theorem states that the LHC would be able to detect at least one of the NMSSM Higgs 
bosons in conventional modes provided Higgs-to-Higgs decays are not important. However, 
Higgs-to-Higgs decays can dominate over other decays, especially for a light Higgs boson, 
because of the small size of fermionic-Higgs Yukawa couplings. This was first established 
in [4] with early subsequent discussions in [2] and [5] followed by studies such as that of [6] 
and [7]. In fact, in [1] it was found that the part of NMSSM parameter space where the 
model has no fine-tuning problem (see later) is such that Higgs-to-Higgs decays must be 
taken into account and do dramatically affect LHC Higgs boson search strategies. In this 
part of parameter space, the lightest CP-even Higgs boson has mass ~ 90 — 100 GeV 
and SM-like gauge and fermion couplings. However, the width for h to decay to a pair 
of the lightest CP-odd Higgs bosons is much larger than the width for h decay to bb, 
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T(h — > aa) S> T(h — > bb), as is easily possible given the very small bbh SM-like Yukawa 
coupling. The a has no tree-level coupling to ZZ and would have escaped LEP searches. In 
addition, the a is very weakly produced in hadron collisions and unlikely to be observable 
at the LHC. The h with ~ 100 GeV would also have escaped detection at LEP even 
though its mass is below 114 GeV. The exact lower bound on with dominant h — > aa 
decays depends upon how the light CP-odd a itself decays. If a — > bb, then the h — > bb 
and h — > bbbb channels must be combined. The lower limit on nih depends somewhat 
on the relative branching ratios for these two channels but is typically of order 110 GeV 
when the latter is dominant [8]. As a result, the non-fine-tuned scenarios (which require 
rrih ~ 100 GeV) only survive LEP limits if m a < 2m},. If a — > t + t~, then the lower bound 
on m/j, is close to 89 GeV [9]. If a — ► 2j, then the only applicable lower bound is that 
where h is assumed to decay hadronically but not necessarily into two jets. This lower 
bound is approximately 82 GeV [10]. Thus, scenarios with ~ 100 GeV and dominant 
h —* aa — » t + t~ or 4j decays are not ruled out by any LEP constraints. 

Although the ~ 100 GeV with dominant h — ► aa decay scenario would a priori 
seem somewhat unlikely, in fact it is preferred. In order to avoid fine-tuning (i.e. to avoid 
tuning GUT scale soft-SUSY-breaking parameters so as to obtain the correct value of mz) 
the superpartners (in particular stops) should not be very heavy. In this case, the model 
automatically predicts that the lightest CP-even Higgs boson should have m^ ~ 100 GeV 
for tan/3 > 10, decreasing slowly to ~ 90 GeV for low tan/3 [11]. In other words, if 
the top squarks are light then the tree-level mass (which is predicted to be ~ mz for 
larger tan/3 values) is only modestly corrected by the top+stop loops. Further, in the 
NMSSM there is a U(1)ji symmetry that, if preserved at low scales, would lead to the 
lightest CP-odd scalar being massless. If this symmetry is present at the GUT scale, small 
soft-SUSY-breaking terms that violate the U(1)r symmetry are automatically generated 
at low scales by renormalization group evolution and one most naturally obtains a light 
a [12]. In fact, as we have already noted, if m a > 2m\ ) so that h — > aa — > 46, then 
LEP excludes < 110 GeV. Thus, for the non-fine-tuned ~ 100 GeV case, we are 
essentially uniquely led to the h — > aa — > t + t~t + t~ or 4j scenarios. Further, the study 
of [11] finds that m a > 2m T (but m a < 2m^) is possible without any tuning of the relevant 
GUT-scale soft-SUSY-breaking parameters whereas m a < 2m T can require some tuning 
(see also the discussion below). This leads us to a preference for the h — > aa — > t + t~t + t~ 
scenario. Finally, it is important to note that ~ 90 — 100 GeV has two additional 
attractive features: (a) the predicted value of BR(h — * aa) > 0.75 leads to a reduced 
value for BR(h — * bb) that is typically such as to explain the 2.3a event excess in the 
e + e~ — > Zbb final state for ~ 98 GeV [1]; (b) such m,h provides the best description of 
precision electroweak data. Thus, the scenario envisioned has excellent phenomenological 
and theoretical properties. 

Our strategy will be to utilize the central exclusive production (CEP) process, pp — > 
p + h + p with h —> aa, which is illustrated in Figure |l|. The outgoing protons carry 
almost all of the incoming beam energy and it is understood that they should be detected 
in purpose built low-angle detectors [13, 14]. These high-resolution sub-detectors would 
permit the four-momentum of the central system to be determined accurately and it is 
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Figure 1: Central exclusive NMSSM Higgs production at the LHC. 



our purpose to exploit this capability to eliminate backgrounds and extract the masses of 
both the h and a on an event-by-event basis. We note that the possibility to use CEP to 
investigate new physics (especially super symmetric new physics) has been much studied 
over the past few years, see [15-18]. 

The pp — > p + 4t + p event rate is proportional to 



and to a m^-dependent luminosity function, which characterises the flux of fusing gluons. 
For those points with minimal electroweak fine-tuning (F < 15 in the notation of [1, ll] 1 ) 
there is a factor ~ 3 variation in the value of r e ff . If one also imposes that the "light-a- 
fine-tuning" parameter G (defined in [12]) obeys G < 25 in order to achieve m a < 2mb 
and BR(h — ► aa) > 0.7 (as required to escape LEP limits) 2 , then the variation of r e g is 
reduced to a factor ~ 2, with the precise amount of variation depending somewhat on m^, 
tan/3 and m a . The dependence of T e Q upon m,h is presented in Figure ^, for two different 
values of tan/3. The points correspond to points in parameter space which have F < 15 
and G < 25 and they are coded according to m a . Scatter plots at these two values of tan /3 
cover the most natural region of NMSSM parameter space, since for larger value of tan (3 it 
is not possible to keep both F and G sufficiently small. We note that m a < 2m T generally 
leads to G > 50, for all values of tan /3, and thus the four-tau decay channel is always most 
natural. The second part of the production cross-section is the gluon luminosity function 
and it falls by a factor ~ 2 in going from = 90 GeV to = 110 GeV [15]. Including 
the explicit l/m| dependence of ( |1.1| ), the production cross-section therefore falls by a 
factor ~ 3 as increases from 90 GeV to 110 GeV at fixed T e g. 

In the present paper we focus our attention on a scenario for which the scalar Higgs, 
h, has mass 92.9 GeV, and the pseudo-scalar Higgs, a, has mass 9.73 GeV. This scenario 
can be achieved with F < 15 and G < 25 for either tan/3 = 3 or tan/3 = 10. For our 
study, we employ a value of r e fj = 0.50 x 10 -4 GeV 2 , i.e. the lowest value that can be 
achieved according to the scatter plots of Figure ^. The value of chosen is in the middle 
of the preferred range of values appearing in Figure ^, implying a mid-range value for the 
gluon luminosity. According to the scatter plots of Figure ^, there is a trend whereby 

1 F < 15 corresponds to fine-tuning of GUT-scale parameters no worse than 6.5%. 

2 G < 25 corresponds to no worse than 4% tuning of the A K and A\ parameters of the NMSSM. 



E e ff = T{h -» gg) 



BR(h -> aa)BR(a -> r+r") 2 
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Figure 2: r c ff in units of 10~ 4 GeV 2 versus mh for tan/3 = 3 (top) and tan (3 = 10 (bottom). Point 
coding: (red) diamonds = 2m r < m a < 7.5 GeV; (green) squares = 7.5 GeV < m a < 8.8 GeV; 
(black) pluses = 8.8 GeV < m a < 2mb. 



larger is associated with larger r e fj, so that this particular choice of parameters gives 
a net production rate for the 4r final state that is not far from minimal. For example, 
for tti/j = 93 GeV one finds in Figure [2] values of r c fj as high as 0.9 x 10 -4 GeV 2 for 
m a > 8.8 GeV. The particular set of m/,, m a and r e fj values we employ in this study- 
corresponds to Scenario 1 of [7] which used the version of NMHDECAY [20] available at 
that time. 3 The h — > aa decay occurs with a branching ratio of 92% and each a decays to 

3 Small improvements in NMHDECAY over the intervening years imply that the output mj,, m a and 
r e ff values used here are slightly changed for the same input parameters used in [7]. 
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r + r~ with a branching ratio of 81%. This implies that h detection in conventional modes 
such as the 77 final state or the r + r _ final state will be highly problematical and a new 
detection mode such as the one we study here should be employed. 

The remainder of this paper is organized as follows. In the next section we briefly 
review the proposal to install forward detectors 420 m either side of the interaction point 
at ATLAS and/or CMS. In Section ^, we discuss the simulation of the signal and relevant 
backgrounds. In Section ||, our attention turns to the task of extracting the signal from 
the variety of backgrounds. We propose an analysis which is heavily track based in order 
to minimize the effect of pile- up. We then include an estimation of the combinatorial (or 
overlap) background caused by pile-up of multiple pp interactions. We take account of 
detector efficiencies, discuss different triggering options and assess the significance of the 
Higgs signal for various possible machine luminosities and detector parameters. In Section 
II we present our conclusions. 



2. The FP420 project 

FP420 is an R&D project with the aim of installing forward proton detectors in the high 
dispersion region 420 m either side of the interaction point at ATLAS and /or CMS [14] . The 
detectors would measure protons from CEP, which have lost energy during the interaction 
and are bent out of the beam by the LHC magnets. There are three issues regarding FP420 
which are relevant to this analysis. Firstly, if the protons are tagged, the longitudinal 
momentum loss, £, of the protons can be measured to an accuracy of 1 — 3% [18]. This 
translates, in this paper, into a measurement of the scalar Higgs mass accurate to 2.1 GeV 
per event. Secondly, the acceptance of FP420 for a central mass of 93 GeV is 28% for 
detectors 5 mm from the beam and 38% for detectors 3 mm from the beam. 4 In this 
analysis we choose the detectors to be 3 mm from the beam, which corresponds to an 
approximate acceptance at IP1 (ATLAS) given by 

0.002 < 6 < 0.018 , 
0.0015 < £ 2 < 0.014 

where the asymmetry arises from the LHC beam optics. Finally, the forward detectors 
will be capable of measuring the time-of-fhght (TOF) of each proton from the interaction 
point to an accuracy of 10 ps. The difference in the TOF measurements translates into 
a measurement of the primary vertex position to an accuracy of 2.1 mm, provided the 
reference timing has negligible jitter. Comparing the vertex implied by FP420 to the 
vertex measured in the central detector is very useful for rejecting backgrounds; a feature 
we exploit in Section S In Section 4.6, we present results also for the case pertaining 



should the 10 ps timing be improved to 2 ps. 

We note that there is also the possibility to use proton detectors installed at 220 m 
from the IP [21,22]. Such detectors are potentially very valuable since they would be close 



4 This is in the context of a proton beam for which abeam = 0.25 mm. There is also a constant dead 
region of 0.5 mm to consider at the edge of the proton detectors. 
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enough to the interaction point to be included in the level one trigger at ATLAS and CMS, 
which have latencies of 2.5 fis and 3.5 fj,s respectively. Unfortunately, they do not add much 
to this analysis due mainly to the low Higgs masses of interest. At CMS, the acceptance 
for asymmetric tagged protons, i.e. one proton measured at 220 m and one at 420 m, 
would be at best 5% for a central mass of 90 GeV [23]. At ATLAS, the corresponding 
acceptance for asymmetric events could be as large as 20% if the active edges of the 220 m 
(400 m) detectors are placed at 2 mm (3 mm) from the beam. However, in this case it 
has recently been appreciated that there would be an attenuation of the acceptance for 
symmetric events (i.e. those which use only the 400 m detectors) due to protons that 
would have been detected at 400 m hitting dead material at the edge of the detectors 
at 220 m. The net effect is that the combined acceptance at ATLAS would increase by 
just 5% after including detectors at 220 m [23]. Furthermore, the OLAP background is 
harder to estimate and likely to be bigger for asymmetric events because of the possible 
acceptance of non-diffractively produced protons in the forward detectors [24]. Finally, 
the mass resolution of the central system would be greater than the 2.1 GeV which one 
obtains using 400 m detectors alone. This smears the signal out over a larger region and 
reduces the significance. For these reasons, the impact of 220 m detectors is likely to be 
very modest and we do not consider them in the remainder of this paper. 

3. Simulating the signal and backgrounds 

We have used the ExHuME Monte Carlo event generator (vl.3.4) [25], modified in order 
to simulate NMSSM scalar Higgs production and the decay h — ► aa. ExHuME faithfully 
implements the theoretical approach to CEP detailed in [15, 19]. 5 The a — > r + r~ decays 
are enforced using the PYTHIA event generator [27] and PYTHIA is also used for the 
subsequent r decays. For triggering purposes, we consider only those events which contain 
at least one muon in the final state. 6 

As the pseudoscalar Higgses are heavily boosted, the signal consists of two back-to- 
back pairs of taus and we therefore focus our attention on those backgrounds which have 
similar topology, i.e. those which look like two jets in the central detector. The central 
exclusive backgrounds are generated using ExHuME; we simulate pp — > p + gg + p and 
pp — > p + bb + p. We do not simulate central exclusive light quark production because these 
backgrounds are suppressed relative to bb by a factor of m^/mf [19]. Hence, if we can show 
that the bb background is small, then we show that all di-quark backgrounds are small. 

A comment on the theoretical uncertainties associated with our predictions for CEP 
(signal and background) is in order, not least because our signal cross-section is small. 
Studies in [17,29] have considered these uncertainties in detail and they indicate a factor 2 
to 3 uncertainty in the total rate. The main uncertainties arise from the lack of a detailed 
understanding of the luminosity of the fusing gluons (which is computed using unintegrated 
gluon densities) and from the difficulty to compute the gap survival factor. The theory has 

including an overall gap survival factor of 3% [26] . 

6 For both signal and background processes that are simulated using ExHuME we use the CTEQ6M 
parton distribution functions [28]. 
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however met with some recent quantitative success: the measurement by CDF of central 
exclusive dijet production [30] is in accord with the predictions of ExHuME. Also, it has 
recently been emphasised that data collected in the early period of LHC running can be 
expected to shed further light on the theory [31]. 

We also consider central particle production arising as a result of 'pomeron fusion' 
(which we usually refer to as 'DPE'). In such events, the leading central jets are accom- 
panied by other particle production. We use POMWIG [32] to simulate this source of 
background. POMWIG is able to produce events of the type pp — > p + jjX + p and the 
version we use [33] incorporates the latest diffractive parton distribution functions from 
the HI experiment at HERA [34]. Due to a limitation in creating a large enough sample of 
inclusive dijet events, we concentrate on bbX production. The final jjX background is ob- 
tained by re-scaling the number of bbX events by the ratio of the jjX to bbX cross-sections. 
Since our analysis is more efficient at eliminating light quark and gluon jets, this procedure 
should lead us to a conservative estimate of the pomeron fusion background. We discuss 
the re-scaling in more detail in Section 4.3. Of course the HI parton densities come with 
some uncertainty (see [34]) and this translates into a modest uncertainty in our prediction 
of the background from this source provided we assume that the diffractive partons are 
universal (i.e. that soft re-scattering effects factorize into the gap survival factor). 

In addition to these direct backgrounds, at sufficiently high luminosity it becomes 
necessary to consider the background from what we shall call overlap processes (OLAP). 
Specifically, we consider the possibility of a threefold coincidence of two single diffractive 
pp — > p+ X events with a generic pp — > X inelastic process. This source of background may 
be important whenever there is a large number of pp interactions in a bunch crossing. The 
inclusive QCD events pp — > X are generated using PYTHIA, with the 'ATLAS tune' to 
Tevatron data [35,36]. For practical purposes, we actually generate pp — > 66 + X inclusive 
scattering and re-scale the final results to match the total 2 — > 2 scattering cross-section; 
for details see Section 4. 5. 7 The forward protons (from single diffraction) are then added 
to the event using the prescription presented in [24], which also allows us to estimate the 
probability of the threefold coincidence as a function of instantaneous luminosity. The two 
protons detected by the 420 m detectors do not originate from the same vertex as the 
primary scatter which produces the muon and this can be exploited to reduce the OLAP 
background. In particular, the pp scattering vertices at ATLAS should be distributed in z 
by a Gaussian of width 4.45 cm and in time by a Gaussian of width 0.18 ns [37]. According 
to the results presented in [24], a rejection factor of 18 (15) should be obtained at low 
(high) luminosity running after demanding that the vertex reconstructed from the forward 
protons' time-of-flight is required to fall within ±4.2 mm of the primary vertex. 

We do not consider OLAP backgrounds from twofold coincidences. It was shown in [24] 
that the largest twofold background — a coincidence between a soft central diffractive 
scattering and an inclusive QCD scattering — was at least a factor of ~ 5 smaller than the 
threefold coincidences. This background would be additionally rejected in this analysis as 
a result of the charged track cut, which is introduced in the next section. 

7 We generate 2 — > 2 scattering processes with prmin = 4 GeV which is sufficiently small given that we 
also require a muon with pr > 6 GeV. 
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Finally, we consider the pure QED backgrounds: pp —* p + t + t~1 + 1~ + p (where I is 
any charged lepton) which is simulated using MADGRAPH [38]. We ignore the effect of 
the proton electromagnetic form factor {i.e. we assume pointlike protons) in the simulation 
but have checked that its effect is not very important in the kinematic regime we consider. 
The final state tau particles are again decayed using the PYTHIA event generator. Note 
that we only consider QED backgrounds which include at least one tau pair: non-tau 
backgrounds can always be eliminated by requiring there to be some missing momentum 



in the central system {i.e. using cuts on the fi variables which we discuss in Section 4.2). 

To improve generator efficiency and restrict ourselves to a broad region of interest, 
we impose an initial condition that the mass of the central system lies in the range 
70 GeV < M < 110 GeV. All final state particle four- momenta are smeared according 
to the relevant detector component resolution [39] with the outgoing proton momenta 
smeared by the amount given in [18]. When appropriate, we also simulate the effects of 
pile-up by superimposing additional inelastic pp collisions simulated using PYTHIA on top 
of both signal and background events. 

4. Extracting the signal 

Our analysis is specifically designed to minimize the impact of extra particles in the cen- 
tral detector which arise from multiple proton-proton interactions in each bunch crossing 
(pile-up). In particular, we propose to utilize mainly tracking information rather than 
information from the calorimeters. Good vertex resolution for charged tracks means that 
the pile-up events rarely contaminate the signal with spurious charged tracks (we quan- 
tify this later) and thus we can select events according to the number of charged tracks. 
We shall select events containing either 4 or 6 charged tracks and at least one muon with 



Pt > 6 GeV. We comment on the effect of a higher muon pr threshold in Section 4.4. In 
this way we collect all signal events in which one tau decays to a muon whilst of the three 
remaining taus at most one of them is allowed to decay to three charged particles. Without 
a transverse momentum requirement for the muon, we would expect 47% of the signal to 
have one or more muons in the final state. With the muon px threshold imposed, we find 
that 25% of the signal is available for further analysis. 

4.1 Event selection 

Before any cuts, the signal is predicted to have a cross-section of 4.8 fb after accounting for 
the branching ratios for the decay to the 4r final state. However for triggering purposes we 
are forced to start from a sample of events (signal and background) containing at least one 
muon with px > 6 GeV. We also select only those events in which the forward protons fall 
within the acceptance of the 420 m detectors. The cross-sections after only these cuts are 
presented in the first line of Table [l]. In the case of the overlap background (OLAP), the 
cross-section is that computed assuming a luminosity of 10 34 cm -2 s _1 and includes the 
factor of 15 suppression which we discussed in the previous section and which represents 
the background rejection from time-of-flight vertexing. We remind the reader that the 
OLAP background is luminosity dependent and its cross-section scales approximately as 



-8- 





CEP 


r>PF 


DT AP 


QED 


Tut 


it 

11 


hh 
uu 


99 


hh 
uu 


hh 
uu 


At 


Li 1 Lit 


"TCP 


0.442 


25.14 


1.51x 10 3 


1 9Qx10 3 


1 74 x 1 6 


01 4 


467 


AT .= 4 or 6 


0.226 


1.59 


28.84 


1 58 x 10 2 


1.44 x 10 4 


0.003 


0.056 


Q& = o 


0.198 


0.207 


3.77 


18.69 


1.29xl0 3 


5xl0" 4 


0.010 


Topology 


0.143 


0.036 


0.432 


0.209 


1.84 




<0.001 


p^,, isolation 


0.083 


0.001 


0.008 


0.003 


0.082 








0.071 


5xl0~ 4 


0.004 


0.002 


0.007 






m a > 2m T 


0.066 


2xl0" 4 


0.001 


0.001 


0.005 







Table 1: The table of cross-sections for the signal and backgrounds. All cross-sections are in 
femtobarns and the cuts are explained in the text. The overlap (OLAP) background is computed 
at a luminosity of 10 34 cm~ 2 s _1 . 
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Figure 3: (a) The expected number of charged tracks reconstructed by the ATLAS inner detector 
with px > 6 GeV and \rj\ < 2.5. (b) The total electric charge in each hemisphere (as defined by the 
trigger muon direction). 



£ 2 ; it is therefore negligible at low luminosities. In this section, we ignore the possibility 
that pile-up can introduce additional charged tracks but we shall return to justify this in 
due course. 

The second line of the table illustrates the effect of demanding only 4 or 6 charged 
tracks. Specifically, for each charged particle, we evaluate the probability that the particle 
is reconstructed as a track by the ATLAS inner detector (ID) and remove particles from the 
event based on this probability. We also only consider those particles with pt > 500 MeV 
and pseudo-rapidity satisfying \ rj\ < 2.5, which is the approximate acceptance of the ATLAS 
ID. All track kinematics are smeared with detector resolution and each track assigned a 
vertex using the ID vertex resolution. Each track vertex is then required to lie within 
2.5 mm of the vertex defined by the muon. This last requirement is important in minimizing 
the impact of charged tracks from pile-up events. The number of charged particles within 
the vertex window, for zero pile-up events, is shown in Figure |3|(a). 
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A smaller vertex window is not desirable because as the vertex window decreases there 
is an increasing possibility that, after vertex smearing, charged particles are assigned a 
vertex outside of the window and would be lost to the analysis thereby reducing the signal. 
Further, because the background events typically have greater than 6 charged particles in 
the primary event, the background can even increase as this window decreases. Of course, 
taking too large a window allows pile-up to contaminate the signal 8 . We choose 2.5 mm 



as a compromise and show the effect of changing the window size in Section 4.5. 

The event is then divided into hemispheres, defined by the axis of the highest-pr 
muon. Figure ^(b) shows the sum of the electric charges in each hemisphere for the signal 
and background events. We assume 100% efficiency in assigning the charge of each track, 
which is a reasonable assumption for reconstructed low transverse momentum tracks as 
the curvature of the track due to the magnetic field will be well measured. The third line 
of Table p] follows from the requirement that there is no net charge in each hemisphere. 

The next cut, the topology cut, requires a little explanation. The first step is to define 
four 'tau objects'. In 4- track events this is trivial but for 6-track events we look to merge 
3 tracks into a single object. The algorithm works as follows: 

1. The nearest two tracks are found using the measure 



r 



ij = (TH-Tl3) 2 + (<f>i-<f>j) 2 (4- 1 ) 



and merged into a single object with a pseudo-rapidity and azimuth given by the pt 
weighted average of the pair, i.e. 

V = (4-2) 

i 

and similarly for <f>. 

2. The nearest track to this pair is found and a tau object defined by clustering the 
3 tracks together. The direction of the corresponding tau object is given by the pr 
weighted average of the pseudo-rapidity and azimuth of the 3 tracks given by ( |4.2j ) . 
The pt of the object is just the sum of the individual track transverse momenta. 

Figure ^a) shows the distance in 7] x <f> space of the three tracks from the newly defined 
tau object and events are only retained if each of the three tracks lies within AR < 0.2 
of the tau object direction. For the surviving events we pair up the four tau objects. 
Figure 0(b) shows the distance in rj x 4> space of each tau object from its nearest neighbor. 
We exploit the fact that the signal has two pairs of roughly collinear tracks by cutting on 
the separation between the tracks which form each pair, i.e. AR < 1. 

The pairing up of the tau objects corresponds to identifying the collinear decay prod- 
ucts of each of the two pseudo-scalars. We refer to each pair as a pseudo-scalar object 



3 It contaminates the background too. 
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Figure 4: (a) The distance in 77 x </> of the three clustered tracks (in a 6 track event) that define 
a single tau object from the weighted average of the tracks, (b) The distance in r\ X <fi of each tau 
object from its nearest neighbor. 



which is constructed using the pt weighted average (77, 4>) for each pair of tau objects. The 
PT of each pseudo-scalar object is just the sum of the pr's of the constituent tau objects. 
Figure |5|(a) shows the azimuthal separation of each pseudo-scalar object, for the signal and 
background events, and we make a cut on A<p > 2.8 to enforce that they are back-to-back. 
Finally, we make the cut 

Ay = i |log(£i/&) - (m +m)\< o.i (4.3) 

where £j is the fractional momentum lost by incoming proton i (as determined by the 420 m 
detectors). This cut exploits the fact that the forward detectors are able to determine the 
boost which will put the central system into the zero momentum frame and that in this 
frame the primary a particles (and hence their approximately collinear decay products) 
have rapidities that are almost equal in magnitude (but opposite in sign). Figure ||(b) 
shows the Ay variable for the signal and background events. The net result of these 
topology cuts is presented in the 4th line in the table. 

We have already insisted upon a single muon with px > 6 GeV. This value is chosen 
since it is the smallest level 1 muon trigger threshold foreseen at ATLAS. Figure ^(a) shows 
the pt spectrum of the trigger muon after the charged track multiplicity cut: it suggests 
we should use a higher cut of px > 10 GeV. In addition to this, we require that the muon 
be isolated, i.e. that there is no hadronic activity near to the muon. Figure ^(b) shows 
the average transverse energy deposited within a distance AR of the muon, after the muon 
Pt > 10 GeV cut. We require that the transverse energy in the calorimeter lying within 
AR < 0.3 of the muon be less than 2 GeV (excluding the fi). This allows for the elimination 
of backgrounds in which the muon is produced as a result of in-jet particle decays, i.e. in 



particular the gg backgrounds. We examine the isolation criterion further in Section 4J5 
since it is sensitive to pile- up. Note that this is the only cut we use which makes use of the 
calorimeter. 
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Figure 5: (a) The difference in azimuth between the reconstructed pseudo-scaiar objects, (b) The 
Ay distribution. 
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Figure 6: (a) The pr of the trigger muon. (b) The average transverse energy deposited in the 
calorimeter within a distance, AR, in 77 x (f> space of the trigger muon. 



We next cut on pj, , which is the transverse momentum of the highest pt track that is 
not the trigger muon. Figure [?] shows the value of p^ 1 for the signal and background events 
before the topology cuts. We choose to present this distribution at this point because the 
statistics are large for each sample. However, in the table the cut is applied after the 
topology cut. The background particles have mainly low px tracks originating from the 
hadronization of jets. The signal tracks however, originate from tau decays and can have 
much larger momentum. We make the cut ph^ > 6 GeV in line 6 of the table. 

4.2 Reconstruction of the Higgs masses 

Since the a bosons are much lighter than the h, it follows that they are highly boosted and 
that their decay products will travel along the direction of the originating a. This means 
that the four-momentum of the primary a is, to a fair approximation, directly proportional 
to the visible four-momentum of the charged particles into which it decays, i.e. 

pT S = fiPa,i (4-4) 
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Figure 7: The pr of the highest momentum track that is not the trigger muon 



where p\ %s is the visible four-momentum. Since we are still restricting our analysis to 
charged tracks, p\ %s is missing the neutrino momentum and the momentum of any neutral 
particles. p a ^ is the four-momentum of the corresponding pseudoscalar a and 1 — is 
the fraction of the a momentum carried away by neutrals. We therefore have only two 
unknowns, f\ and fi and these can be determined using information from the forward 
proton detectors since we know that 



Pa,X + Pa,2 = Ph 



(4.5) 



and ph can be measured. In fact the situation is over-constrained since we have only 
two unknowns and three equations. Although the transverse momentum of the Higgs 
can in principle be measured using the forward detectors it will typically be rather small. 
Assuming it to be zero leaves us with the three equations: 



(pT 



f: 



)x,y {P™ S )x,y 



f2 







and 



(P 



+ 



(P: 



ms^ 



(4.6) 



(4.7) 



fx h 

where x and y label the directions transverse to the beam axis and the 1 — & are the longi- 
tudinal momenta of the outgoing protons expressed as fractions of the incoming momenta. 
Solving Q and Q gives 



fx 

h 



(p\ 



(pf s UP?%,y 



(Pi 



VIS\ 

2 )x,y 



fx ■ 



(4.8) 
(4.9) 



(P?%,y 

Equations ( f4.8[ ) and (^) provide two solutions depending on whether we solved using the 
(x, z) or (y, z) pair of equations. We therefore have 4 determinations of the a mass per 



event. - 



9 The (x, y) pair of equations could also be used if one could make a good measurement of the outgoing 
proton transverse momenta. 
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(a) (b) 

Figure 8: (a) The ratio of the reconstructed scalar Higgs mass to the mass measured by FP420 for 
the signal events only, (b) The reconstructed a mass for the signal events. The broad distribution 
is due to the breakdown of the collinearity approximation and detector effects are minimal. 



Using this method, we are able to reconstruct the masses of both the scalar Higgs 
and the pseudo-scalar Higgs. Figure |8|(a) shows the ratio of the reconstructed scalar Higgs 
boson mass to the mass measured in the forward detectors for the signal. The distribution 
is broad, mainly as a result of the collinearity approximation and the missing momentum 
carried by neutral particles (this is a bigger effect than that due to detector smearing). In 
the final analysis, we do not need to apply a cut on this variable because we have already 
adequately reduced the background. In Figure |§|(b), we show the pseudo-scalar mass 
distribution using only signal events. The final cut in Table [l] corresponds to restricting 
the mean of each of the two reconstructed a masses to the range m a > 2m T . 

Finally, it should be noted that we have not yet applied a stringent mass window re- 
quirement. The background was generated with a central mass in the range 
70 GeV < M < 110 GeV and the mass distribution as measured by FP420 (with de- 
tector smearing) is shown in Figure ^ after the application of the topology requirement. It 
is clear that in estimating the significance we should further limit ourselves to a smaller 
region of interest around the signal. 

4.3 Estimating non-66 backgrounds 

In Table [l] the OLAP background quoted is obtained by superimposing an inclusive pp — > 
bb + X event with two single diffr active events pp — > p + X. As stated, we obtain the 
full OLAP background after re-scaling the OLAP background from this source by a factor 
a(pp — > jj + X)/a(pp — ► bb + X). PYTHIA is used to generate the events and the cross- 
sections are obtained after insisting that there be at least one muon with pr > 6 GeV, and 
after the proton acceptance and the charged track cuts: we find that a re-scaling of the bb 
OLAP background by a factor 5.0 is necessary. 

The re-scaling is determined after the charged track cut since many of the muons 
produced in the jj + X sample arise from processes involving gluon jets (this is certainly so 
in the case of gb — > gb intrinsic b contributions) in which case there are typically more tracks 
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Figure 9: The mass distribution of events after the topology requirement for signal and selected 
background. 

than in a typical bb event and hence the track cut should be more effective at reducing this 
source of background. The results quoted in the following sections include this re-scaling. 

A similar re-scaling is performed, this time using POMWIG, in order to convert the 
DPE bb background to the full jj background and we find that a correction factor of 1.86 is 
required in this case. This number is smaller than for the inclusive case since the pomeron 
is gluon dominated: In the inclusive case, intrinsic c and b quark contributions actually 
dominate the production of those dijet events which pass our muon pj- and charged track 
cuts, via the process gq — ► gq. 

One may also worry about OLAP coming from processes other than QCD dijet pro- 
duction. For example, the inclusive QED production of four taus in conjunction with two 
single diffractive events. However, this background is negligible since we know that the ex- 
clusive QED background is negligible, as demonstrated in Table 1 (and this was computed 
without any form factor suppression). The fact that the photons can couple to quarks 
within the protons is not able to compensate for the suppression factor of over two orders 
of magnitude arising from the requirement that there also be two single diffractive events 
which produce a fake vertex coincident with the vertex from the muon. 

4.4 Triggering 

In this section, we discuss the trigger options and trigger efficiencies for this analysis. The 
signal from FP420 will not arrive at the central trigger processor in time for the level 1 
decision. At level 2 and above however, we can make use of the FP420 and the muon vertex 
information to substantially reduce the rate; requiring two in-time proton hits results in a 
reduction factor of 140 at high luminosity and approximately 20000 at low luminosity [24] . 
We therefore focus on the possible ways to retain the events at level 1. 

As we already stated, our strategy is to trigger on a muon arising from the decay of 
one of the taus. Before any analysis cuts, we imposed a minimum muon pr requirement of 
6 GeV. This is the lowest trigger threshold foreseen at ATLAS and it would be pre-scaled 
at high luminosity in order to reduce the rate and may even be pre-scaled at low luminosity. 
However, in our analysis we make a further cut on the muon; requiring pr > 10 GeV in 
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order to improve the significance of any observation. In our subsequent analysis, we will 
therefore employ what we call the MU10 trigger, i.e. a muon trigger threshold of 10 GeV. 
For muons with pp > 10 GeV, the trigger efficiency is 80% [40,41]; in this way we avoid 
the need to pre-scale without losing too much signal. It should be noted that we still need 
to employ the pt > 10 GeV analysis cut since the MU10 trigger does not directly measure 
the muon pp, which results in the acceptance of some muons with pp below 10 GeV. The 
bottom line is that, for the MU10 trigger, we multiply the event rates in Table |l] by 0.8 to 
reflect the ATLAS trigger efficiency 

The fact that the ATLAS level 1 muon trigger has acceptance below the nominal 
threshold allows us a second possibility to reduce rates: we could use a MU15 trigger. 
Such a trigger would have an efficiency of 80% for muons with pp > 15 GeV, but would 
retain some efficiency for pp values in the region 10 — 15 GeV. In [41], results are presented 
for the LI muon trigger efficiencies for 11 GeV and 20 GeV thresholds and from these 
results it is apparent that the efficiency for triggering on 10 GeV muons given a 15 GeV 
threshold will be between 20% and 50%. We choose to model the efficiency of a MU15 
trigger with a linear turn-on, with 20% efficiency at 10 GeV and 80% efficiency at 15 GeV. 
Above 15 GeV we assume the canonical 80% efficiency. This linear turn-on will probably 
underestimate the true efficiency because the real trigger should be more efficient in the 
region near to the threshold (and because the trigger could in fact be closer to 50% efficient 
for 10 GeV muons). The cut on muon pp at analysis level can therefore remain at 10 GeV 
even with a MU15 trigger. We note that any selection bias due to a different rj — <f> 
dependence of the signal and background is expected to be small since they have rather 
similar behaviour in these variables. In any case, this small effect would not affect our 
estimate of the significance provided one can determine the background by extrapolating 
from outside the signal region. 

4.5 Contamination from pile-up 

So far we have assumed that the contamination of both signal and background events from 
pile-up is negligible. In this section, we quantify the effect of pile-up on our estimates by 
adding events generated by PYTHIA to our signal events. The first effect is that charged 
tracks from the pile-up events do cause signal events to fail the charged track cut when the 
pile-up vertex is near to the primary vertex. Figure |i~0|(a) shows the fractional reduction 
in signal for various Az choices due to the extra charged tracks from pile-up as a function 
of the number of pile-up events overlaid with the signal event. The cross-section is relative 
to that obtained without pile- up for the given Az and is evaluated after the p T ^ cut. The 
effect of changing the vertex window size is also shown. Note that, as discussed in Section 



4.1, we choose to use \Az\ < 2.5 mm in order to avoid the increase in background, which 
arises as \Az\ decreases, and also to avoid having too many signal tracks start to fall outside 
the window. 

The other effect of pile-up is that the additional energy in the calorimeters can cause 



events to fail the muon isolation cut. Figure 10(b) shows the reduction in signal events due 



to the increased energy in the calorimeter as a function of the number of pile- up events. 
It is clear that requiring Et < 2 GeV in the calorimeter within AR < 0.3 of the muon is 
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(a) (b) 

Figure 10: (a) The fractional reduction in the signal cross-section (evaluated after the pj^ cut) as 
a function of the number of pile-up events for three different vertex window sizes, (b) The fractional 
reduction in the signal cross-section as a function of the number of pile-up events for two different 
muon isolation criteria. 



only valid at low luminosities. Figure 10(b) also shows the effect of imposing a luminosity 
dependent isolation cut given by 

N 

E T < 2.0 + — (GeV) (4.10) 

where N is the average number of pile-up events in a bunch crossing. This results in an 
isolation requirement of 5.4 GeV at high luminosity, which retains approximately 50% of 
the signal. In what follows we use Eq. fl4.10| ). 



4.6 The significance 

In this section we make predictions for the significance of any observation given three years 
of data acquisition at a fixed luminosity. We take the nominal integrated data acquisition 
rate of 10 fb _1 yr _1 at low luminosity (10 33 cm~ 2 s" 1 ) and 100 fb _1 yr _1 at high luminosity 
(10 34 cm -2 s _1 ). In order to predict final event rates, we must re-scale the results presented 
in Table [l] as follows: 



The DPE and OLAP bb backgrounds are scaled up as explained in Section 4.3 



All event rates are scaled to reflect ATLAS trigger efficiency as explained in Sec- 



tion 4.4 



All event rates are scaled to reflect the effect of pile-up on the charged track multi- 
plicity and the muon isolation cuts, i.e. using the solid black curves in Figures |l0|(a) 
and (b); 

For significance purposes, we are interested in the event rates in a region around the 
Higgs mass peak of rrih ± 5 GeV. We reduce our background rates accordingly. 
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0.08 
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3.3 


0.36 
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0.11 



Table 2: Expected number of signal (S) and background (B) events for the three trigger scenarios 
assuming that the data are collected at a fixed instantaneous luminosity over a three year period. 
We assume the integrated luminosity acquired each year is 10 fb _1 , 50 fb _1 and 100 fb _1 at an 
instantaneous luminosity of lxlO 33 cm~ 2 s _1 , 5xl0 33 cm~ 2 s _1 and 10 xlO 33 cm~ 2 s _1 . 
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Figure 11: (a) The significance for three years of data acquisition at each luminosity, (b) Same 
as (a) but with twice the data. 



Table ||| shows the final number of signal and background events after three years of 
data acquisition at an instantaneous luminosity of 10 33 , 5 x 10 33 and 10 34 cm" 2 s" 1 . We 
present the final event numbers for the two trigger strategies, MU10 and MU15, and also 
for MU10 with reduced overlap backgrounds, which assumes that 2 ps proton time-of-flight 
can be achieved by FP420. 

We estimate the significance, S, of the observation using 

(4.11) 



1 
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e 2 dx 
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n=s+b 



where s is the number of signal events and b is the number of background events. The 
probability of observing s + b events given a Poisson distribution with mean b is equated 
to the probability of S standard deviations in a Gaussian distribution. 

Figure |ll](a) shows the significance of the measurement after three years of data ac- 
quisition at each luminosity. The lower curves correspond to the two different muon pt 
trigger thresholds whilst the upper curve shows what it possible if the forward detector 
timing accuracy can be improved from 10 ps to 2 ps. Figure 0(b) is the same as Fig- 
ure |ll|(a) but assuming that twice the amount of data are available, e.g. as might occur if 
the results of ATLAS were combined with those of CMS. 
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Figure 12: The significance of the measurement for reduced overlap (OLAP) contributions. 



Figure 12 shows the significance of the measurement if the overlap background was 
increased or decreased by a factor of 5, assuming the nominal data acquisition rate and 
10 ps time-of-flight measurements. This plot illustrates the importance of the OLAP 
background at higher luminosities. We observe that a five-fold increase in the OLAP 
background would result in the analysis being possible only at low luminosity. A similar 
plot for the case of twice the data being available shows that the OLAP background would 
have to be more than a factor of 40 larger than our best estimate before the signal would 



become unobservable. We show Figure 12 because we do not believe it is possible to make a 
precise prediction for the OLAP background given current understanding. In particular we 
note that the huge reduction factors for this background arise mainly because we make cuts 
on the charged track multiplicity, and this is heavily dependent on the underlying event. 
In this analysis, we chose to use the 'ATLAS tune' of PYTHIA which fits the Tevatron 
data. However, the extrapolation of these fits to LHC energies is uncertain and different 
models/tunes predict a different amount of underlying event activity. For example, if we 
had chosen to use the HERWIG event generator [42], with JIMMY [43] simulating the 
underlying event, the charged track activity would be increased [35] which would make 
the overlap events more likely to fail the charged track multiplicity cuts thereby reducing 
the OLAP background. Apart from uncertainties arising due to our lack of understanding 
of the OLAP background, our estimates have assumed that all of the data are collected 
at a value of the instantaneous luminosity fixed to the canonical value (10 33 cm~ 2 s _1 or 
10 34 cm _2 s _1 ). In reality, the luminosity decreases from the canonical value with time until 
the next LHC fill. The decay time of the luminosity at the LHC is expected to be 14.9 
hours [44] and the ratio of the average luminosity per fill to the initial luminosity is expected 
to be 0.69 for 12 hour runs and 0.82 for 6 hour runs. Thus the OLAP background could 
easily be up to two times smaller than that calculated in Table [|, even at high luminosity 
running. 

To conclude, we show in Figure |i~3|(a) a typical distribution in the measured mass of the 
light pseudo-scalar a. Since we make 4 mass measurements per event, the 6 signal events 



each contribute 4 entries in the histogram. Figure 12 (b) shows the distribution of masses 
obtained for each event; the integer in each box labels one of the 6 signal events. The figure 
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Figure 13: (a) A typical a mass measurement, (b) The same content as (a) but with the breakdown 
showing the 4 Higgs mass measurements for each of the 6 events, labelled 1 — 6 in the histogram. 



is indicative of what one could obtain with 180 fb _1 of data collected at 3 x 10 33 cm _2 s _1 . 
By considering many pseudo-data sets, we conclude that a typical experiment would yield 
m a = 9.3 ± 2.3 GeV, which is in re-assuringly good agreement with the expected value of 
9.7 GeV. 

5. Conclusions 

Both the Standard Model and Minimal Super symmetric Model have theoretical problems 
with regard to naturalness, hierarchy and fine-tuning, especially after including the LEP 
constraints on the lightest CP-even Higgs boson. Models which avoid these problems are 
typified by unusual decay modes for the lightest of the CP-even Higgs bosons. As a result, 
there is no guarantee that LHC Higgs discovery will be possible using the conventional 
modes explored for the SM and MSSM. A particularly attractive model that avoids the 
theoretical problems of and experimental constraints on the SM and MSSM is the Next- 
to-Minimal Supersymmetric SM. Hierarchy and fine-tuning issues are absent provided the 
soft-SUSY-breaking parameters are chosen in such a way that sparticles are relatively 
light, leading to a light CP-even h with ~ 100 GeV and SM-like couplings to gauge 
bosons and fermions. The h of the NMSSM must have escaped LEP detection by virtue 
of h —* aa —* 4r or 4j, with the former being the most natural result in the context of 
the theory. We have demonstrated that the prospects of utilizing CEP to observe such 
an NMSSM h decaying via h —* aa — > 4r, are good, provided that sufficient integrated 
luminosity is available. The viability of the pp — > p + h + p detection channel is especially 
important given that no other LHC modes have as of yet been shown to be viable and 
given the possibility to extract the mass of the h and a on an event-by-event basis. 

With regard to the need for high instantaneous luminosity, we have demonstrated that 
the analysis can work even at high instantaneous LHC luminosity, with the optimum being 
about C ~ 5 x 10 33 cm~ 2 sec" 1 . This is because of our ability to bring pile-up under 
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control, which in turn is a consequence of the fact that our analysis makes very little use 
of the hadronic and electromagnetic calorimeters and instead relies more on tracking and 
the muon detectors. 

Although we have not explored the matter in this paper, it seems quite likely that 
such a track-based analysis could also be utilized in other CEP processes, for example 
pp — > p + h + p with h — ► t + t~. A heavily track-based approach could also provide an 
interesting analysis route for more mainstream LHC processes in order to help diminish the 
influence of pile- up, e.g. in analyses which aim to exploit the presence of rapidity gaps. An 
obvious example of such a process is Higgs production via weak boson fusion (WBF). WBF 
is traditionally distinguished from Higgs production via gluon fusion (and backgrounds) by 
requiring the presence of two highly energetic jets at large absolute rapidity values (recoil- 
ing from the emitted virtual Ws). Since the virtual Ws are color singlets, the rapidity 
region between the recoil jets and the centrally produced Higgs decay products should 
have reduced jet activity compared to backgrounds and other Higgs production processes. 
An analysis similar to that presented here could be used to select events according to the 
numbers of tracks they contain in the central region. This approach would be particularly 
interesting for H — > t + t~ (in the SM or otherwise) where one r decays to a muon. The 
muon then defines, very precisely, the location of the primary vertex and one could cut on 
the number of additional tracks pointing back to that vertex. It will also be important to 
explore whether an analysis along these lines might even make viable the detection of the 
NMSSM h — ► aa — ► 4t Higgs decay mode when the h is produced via WBF. Whatever the 
primary Higgs decay mode, but especially in the h —* aa — > 4r case, observations of the 
same Higgs final state via WBF would be very important for confirming any Higgs signal 
seen in CEP. 

Finally, we return to the fact that m a < 2m T remains possible in the NMSSM, although 
this requires a little more parameter tuning. Thus, there is a need to identify techniques 
that might yield a viable signal for the resulting h — > aa — > Aj final state. Similarly we 
should keep in mind that rrih > 110 GeV with large branching ratios for h — > aa and a — > bb 
is not excluded by LEP. It is therefore also important to establish the viability of observing 
this channel in CEP. 
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